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Abstract. The first static spherically symmetric perfect fluid solution with constant density 
was found by Schwarzschild in 1918. Generically, perfect fluid spheres are interesting because 
they are first approximations to any attempt at building a realistic model for a general 
relativistic star. Over the past 90 years a confusing tangle of specific perfect fiuid spheres 
has been discovered, with most of these examples seemingly independent from each other. To 
bring some order to this collection, we develop several new transformation theorems that map 
perfect fluid spheres into perfect fiuid spheres. These transformation theorems sometimes lead 
to unexpected connections between previously known perfect fiuid spheres, sometimes lead 
to new previously unknown perfect fiuid spheres, and in general can be used to develop a 
systematic way of classifying the set of all perfect fiuid spheres. In addition, we develop new 
"solution generating" theorems for the TOV, whereby any given solution can be "deformed" 
to a new solution. Because these TOV-based theorems work directly in terms of the pressure 
profile and density profile it is relatively easy to impose regularity conditions at the centre of 
the fluid sphere. 



1. Introduction 

Perfect fluid spheres, either Newtonian or relativistic, are the first approximations in developing 
realistic stellar models El El- Because of the importance of these models, fully general 
"algorithmic" solutions of the perfect fluid constraint in general relativity have been developed 
over the last several years jU El El 13 IHl- For some of the key highlights of earlier history 
regarding perfect fluid spheres see references ^ IQ., JJ^ , while for a recent analysis providing 
upper and lower bounds on the internal structure of a relativistic perfect fluid sphere one might 
consult [T3| . 

For our current purposes, the central idea is to start solely with spherical symmetry, which 



implies that in orthonormal components the stress energy tensor takes the form: 
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and then use the perfect fluid constraint pr = Pt- This simply makes the radial pressure equal to 
the transverse pressure. By using the Einstein equations, plus spherical symmetry, the equality 
Pr = Pt for the pressures becomes the statement 

Gq§ = Gff = G^^. (2) 

In terms of the metric components, this leads to an ordinary differential equation [ODE], which 
then constrains the spacetime geometry, for any perfect fluid sphere. 

The big recent change has been the discovery and development of "algorithmic" techniques 
that permit one to generate large classes of perfect fluid spheres in a purely mechanical 
way El 13 EI- In this article we will summarize and discuss our recent work on these 
algorithmic ideas IE]) providing several solution-generating theorems of varying levels of 
complexity. We shall then explore the formal properties of these solution-generating theorems, 
and shall then use these theorems to classify some of the previously known exact solutions. 
Moreover, we shall then generate several new previously unknown perfect fluid solutions, and 
finally we shall show how one can rephrase all these theorems directly in terms of the pressure 
profile and density proflie. 



2. Solution generating theorems 

Start with some static spherically symmetric geometry in Schwarzschild (curvature) coordinates 

ds' = -C{rfdt' + ^+r'dn^ (3) 

and assume it represents a perfect fluid sphere. That is, Ggg = Gff = G^^. While G^g = G^^ 
is always fulfilled due to spherical symmetry, setting Gff = Gqq supplies us with an ODE 

[r{rQ)']B' + [2r^C" " 2(rC)']S + 2C = , (4) 

which reduces the freedom to choose the two functions in equation ^ to one. This equation 
is a first order-linear non-homogeneous equation in B{r). Thus — once you have chosen a Cii^) 
— this equation can always be solved for B[r). Solving for B{r) in terms of C,{r) is the basis of 
the analyses in references (and is the basis for Theorem 1 below). On the other hand, 

we can also re-group this same equation as 

2r'^BC" + {r^B' - 2rB)Q' + {rB' - + 2)C = , (5) 

which is a linear homogeneous 2nd order ODE for C('")) which will become the basis for Theorem 
2 below. The question we are going to answer in this section is, how to systematically "deform" 
this geometry while still retaining the perfect fluid property. That is, suppose we start with 
the specific geometry defined by 

ds^ = -C,{rfdt^ + 4T-^+r^d^l' (6) 
Bo[r) 



and assume it represents a perfect fluid sphere. We will show how to "deform" this solution 
by applying four different transformation theorems on {(^q,Bq}, such that the outcome still 
presents a perfect fluid sphere. The outcome of this process will depend on one or more free 
parameters, and so automatically generates an entire family of perfect fluid spheres of which 
the original starting point is only one member. We also try to find the connection between all 
different transformation theorems 

2.1. Four theorems 

The first theorem we present is a variant of a result first explicitly published in reference [H]. 
Theorem 1 Suppose {C,Q{r) , BQ{r)} represents a perfect fluid sphere. Define 



Xo{r) + r C'^{r) J J CoW CoW + r Co(r 

Then for all A, the geometry deflned by holding Co('") fixed and setting 

ds^ = -Ur? dt^ + + r^dQ^ (8) 

Bo(r) + A Ao(r) 

is also a perfect fiuid sphere. That is, the mapping 

Ti(A) : {Co,i3o}^{Co,i?o + AAo(Co)} (9) 

takes perfect fluid spheres into perfect fluid spheres. 

Proof for Theorem 1 Assume that {C,Q{r) , BQ[r)} is a solution for equation Q). Under what 
conditions does {C,Q{r),Bi{r)} also satisfy equation Q? Without loss of generality, we write 

fii(r) = 5o(r-) + A Ao(r). (10) 

Equation @ can now be used to determine Ao(r). That ordinary inhomogeneous first-order 
differential equation in B now simplifies to 

[r(rCo)'] A'o + [2rX^ - 2(rCo)'] Aq = , (11) 

which is an ordinary homogeneous first-order differential equation in Aq. A straightforward 
calculation, including an integration by parts, leads to 

^o(r) = exp I / j^dr] . (12) 



[(rCo)f U Ko)'' 
Adding and subtracting zb2(r^Q^)/(Co('"Co)') to the argument under the integral leads to 

\Co[r) + r Q{r) J [ J Co{r) Co[r) + r Qir) J 

as advertised. 



Comments: A version of Theorem 1 can also be found in j0J. Specifically, after several 
manipulations, changes of notation, and a change of coordinate system, the transformation 
exhibited in equation (16.11) of [H] can be converted into Theorem 1 above. 



Theorem 2 Let {Co^^o} describe a perfect fluid sphere. Define 

r dr 



Zo{r) = a + e [ — 



(14) 



,(r)2 /B^ 

Then for all a and e, the geometry defined by holding BQ{r) fixed and setting 

ds-" = -Coirf Zoirf df + ^^+ r'dQ' (15) 

is also a perfect fluid sphere. That is, the mapping 

T2(<J,e) : {Co,Bo} ^ {Co Zo{Co, Bo), Bo} (16) 
takes perfect fluid spheres into perfect fluid spheres. 

Proof for Theorem 2 The proof of theorem 2 is based on the technique of "reduction in order" . 
Assuming that {C,o{t^) ■, Bo{r)} solves equation ©, write 

Ci(r) = Co(r)^o(r). (17) 

and demand that {C,i{r) , Bo{r)} also solves equation 0. We find 

{r^oB'o + ^r\'oBo - 2rCoSo)^o + {2r\oBo)Z'^ = , (18) 

which is an ordinary homogeneous second-order differential equation, depending only on Zq 
and Zq. (So it can be viewed as a first-order homogeneous order differential equation in Z' , 
which is solvable.) Separating the unknown variable to one side, 

4 = --^-2^ + -. (19) 
Z'q 2Bo Co r ^ ' 

Re-write Zq/Zq = dln{ZQ)/dt^ and integrate twice over both sides of equation (fl9|l . to obtain 

/T dv 
; , (20) 

depending on the old solution {Co('^)i ^o(^)}i and two arbitrary integration constants a and e. 

Having now found the first and second generating theorems makes it possible to define two 
new theorems by composing them. Take a perfect fluid sphere solution {Co,^o}- Applying 
Theorem 1 onto it gives us a new perfect fluid sphere {Qo^Bi}. The new Bi is given in equa- 
tion 0. If we now continuing with applying Theorem 2, again we get a new solution {C,-Bi}, 
where C, now depends on the new Bi. All together we can consider this as a single process, by 
introducing the following theorem: 



Theorem 3 If {Qq^Bq} denotes a perfect fluid sphere, then for all cr, e, and A, the three- 
parameter geometry defined by 

^ _^^(^)2 <) ^ + e / ^ 1 + ^ , , , + r^dn" (21) 




Co(r)2 ^i?o(r) + A Ao(r) J 5o(r) + AAo(r 

is also a perfect fluid sphere, where Aq is 

^„(r)=l,rrP^^\'<^^A^^W■^M^^.<ir\. (22) 



.Co(r) + rC;(r)y "I ] CoW CoM + r CJM 

That is 

T3 = T2oTi : {Co,5o} ^ {Co,5o + A Ao(Co)} ^ {Co ^o(Co, 5o + AAo(Co)), 5o + AAo(Co)} (23) 
takes perfect fluid spheres into perfect fluid spheres. 

Now, instead of starting with theorem 1 we could first apply theorem 2 on {^Oi-^o}- This 
gives us a new perfect fluid sphere {Ci, -Bo}) where Ci = Co ■^o(Co; Bq) is given by equation (PU]) . 
Continuing with theorem 1 leads to {Ci,-B}, where B depends on the new Ci- Again, we can 
consider this as a single process, by introducing the following theorem: 



Theorem 4 If {Co^^o} denotes a perfect fluid sphere, then for all a, e, and A, the three- 
parameter geometry defined by 

,2 _ A r^2 !_,,/" ^ I J^2 , , „2 .r>2 



ds^ = -UrY <(T + e == } dt' + — — — + r^dn^ (24) 

is also a perfect fiuid sphere, where Ao(Ci;'^) is defined as 

^" VCi(r)+rC;(r)j '"^tV Ci(r) Ci{r)+rC[{r) ^7 ' ^ ^ 

depending on Ci = Co ^Oi where as before 

Zo{r)=a + e [—-^^^ (26) 

That is 

T4 = Ti o T2 : {Co, ^o} ^ {Co ^o(Co, So), 5o} ^ {Co ^o(Co, So), So + AAo(Co ^o(Co, So))} 

(27) 

takes perfect fluid spheres into perfect fluid spheres. 

Below we present two tables which describe the inter-relationships of certain examples 
of perfect fluid spheres under T]^, T2, T3, and T4. We define a seed metric as one for 
which Theorem 1 and Theorem 2 both yield output metrics distinct from the input metric: 
"'"1(9) ^ 9 ^ '^2^9)- contrast for non-seed metrics one or the other of these theorems 
is trivial, either T-^{g) = g or T2{g) — g (where g is the metric or a parameterized class 
of metrics). Classifying perfect fluid spheres into seed and non-seed metrics is a very useful 
step. Our theorems are effectively a device to generate new solutions for a perfect fiuid sphere 
without explicitly solving the Einstein equations. (Of course our theorems are set up in such a 
way that one is implicitly solving the Einstein equations.) The terminology used in the tables 
is more fully explained in reference 0. 



Some seed geometries and their descendants. 



Seed 


Theoreml 


Theorem2 


Theorem 3 


Theorem 4 


Minkowski 


Einstein static 


K-0 III 


int Schwarzschild 


Martin 3 


ext Schwarzschild 


Kuch68 II 


Kuch 86 I 


[integral] 


[integral] 


de Sitter 


Tolman IV 


int Schwarzschild 


P7 


[integral] 


Tolman V (A = 0) 


Tolman V 


Tolman VI 


Wyman III 


Wyman Ila 


SI 


Tolman V (n = +1) 


K III 


P4 


Martin 3 


M-W III 


Martin 2 


P5 


P6 


[integral] 


Heint Ila (C=0) 


Heint Ila 


P3 


P8 


[integral] 



Table 1. Seed solutions and their generalizations derived via theorems 1-4. The notation 
"[integral]" denotes a metric so complicated that explicitly writing out the relevant integral is 
so tedious that it does not seem worthwhile. The terminology is more fully explained, and the 
relevant metrics explicitly presented, in reference [7j. 



Some non-seed perfect fluid geometries and their descendants. 



Base 


Theoreml 


Theorem2 


Theorem 3 


Theorem 4 


Tolman IV 
B-L 

Heint Ha 


Tolman IV 
B-L 

Heint Ha 


P7 

Kuchb I b 
P8 


P7 

Kuchb I b 
P8 


[integral] 
[integral] 
[integral] 


Tolman VI 
Kuchl lb 
K HI 


Wyman Ha 
Martin 1 
Martin 3 


Tolman VI 
Kuchl lb 
K-0 HI 


[integral] 

P2 

PI 


[integral] 
Martin 1 
Martin 3 



Table 2. Non-seed solutions and their generalizations. The terminology used in this table is 
more fully explained in reference jjj. 



3. Solution generating theorems for the TOV equation 

The Tolman-Oppenheimer-Volkov [TOV] equation constrains the internal structure of general 
relativistic static perfect fluid spheres [S|. We have developed several "solution generating" 
theorems for the TOV, whereby any given solution can be deformed to a new solution. Because 
these theorems work directly in terms of the physical observables — the pressure profile 
and density profile — it is relatively easy to examine the density and pressure profiles for 
physical reasonableness. This section complements our previous theorems wherein a similar 
"algorithmic" analysis of the general relativistic static perfect fluid sphere was presented 
directly in terms of the spacetime geometry — in the present analysis the pressure and density 
are primary and the spacetime geometry is secondary. Specifically, consider the well-known 
TOV system of equations, whose derivation is now a common textbook exercise jl4( I15j: 

dp{r) [p{r) + p{r)] [m{r) + Anp^r) r^] 

dr r'^[l — 2m{r)/r] ' 

dm{r) 



dr 



A-Kp{r)r\ (29) 



We adopt units where G = c = 1 so that m/r and pr^ are dimensionless. Our basic strategy 
wiU be to assume that somehow we have obtained, or been given, some specific "physically 
reasonable" solution of the TOV equation in terms of the profiles poir) and p{){r). 

Viewed as a differential equation for p{r), with m{r) [and hence p(r)] specified, the TOV 
equation is a specific example of a Ricatti equation, for which the number of useful solution 
techniques is unfortunately rather limited. Using standard results (see the explicit discussion 
in reference [Hj , or for example the general references |16| 117) ) it is relatively simple to present 
the following: 

Theorem PI Let pQ{r) and PQ{r) solve the TOV equation and hold rriQir) = 4tt J pQ{r)r'^dr 
fixed. Define 

, . mo(r) +47rpo(r)r^ 
''^'^ = r2[l-2mo(r)/r] " ^^^^ 

Then the general solution to the TOV equation is p{r) = po(^) + Spir) where 

6pc a/1 - 2mo/r exp S -2 / godr\ 
Sp{r) = J, ( r \ , (31) 



1 + A-Kdpc / — = exp s — 2 / godr > r dr 

Jo i/l-2mo/r I Jo J 

and where 6pc is the shift in the central pressure. 

A second theorem can be obtained by looking for correlated changes in the mass and pressure 
profiles: 

Theorem P2 Let po{r) and po{r) solve the TOV equation and hold go fixed, in the sense that 

mo{r) + 'i7rpo{r)r'^ m(r) + 47rp(r)r'^ 
^^^'^^ " r2[l-2mo(r)/r] " r^[l - 2m{r)/r] ' ^ ' 

Then the general solution to the TOV equation is given by p{r) = po{r) + 5p{r) and 
m(r) = mo{r) + 5m{r) where 

= e.p (2 r» . (33) 



3 [l + rgo]'^ [Jo I + r go 

and 



^P^'^=-4^^ l-2mo/r - ^^^^ 



Here 5pc is the shift in the central density. 



4. Discussion 

Using Schwarzschild coordinates we have developed two fundamental transformation theorems 
that map perfect fiuid spheres into perfect fluid spheres. Moreover, we established two 
additional transformation theorems by composing the first and second generating theorems. 
We then used these transformations as a basis for classifying different types of perfect fluid 
sphere solutions. If we apply these theorems on a known perfect fluid sphere, different solutions 
are often obtained. While some of the solutions are already known in the literature, most of 



them are novel. Therefore, we have devefoped a tool to generate new solutions for a perfect fluid 
sphere, but which does not require us to directly solve the Einstein equations. In principle all 
solutions of the perfect fluid constraint should be obtainable in this manner. In short, perfect 
fluid spheres may "simple" , but they still provide a surprisingly rich mathematical and physical 
structure. 

Furthermore, we have also developed two "physically clean" solution-generating theorems 
for the TOV equation — where by "physically clean" we mean that it is relatively easy to 
understand what happens to the pressure and density proflles, especially in the vicinity of the 
stellar core. 

Finally, we reiterate that the relativistic perfect fluid sphere, despite its venerable history, 
continues to provide interesting surprises. 
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